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TECHNICAL MEMORANDUM NO. 962 



RECENT TTORE ON AIRFOIL THEORY 
By L.' Prandtl 

I should like in the following to report "briefly on 
several papers which have appeared in Gottin^en during the 
last three years. 

I. In the computational treatment of the lifting sur- 
face, progress has hoen made "by starting out - not from a 
bound vortox distribution on the surface with the associat- 
ed trailing vortex sheet - hut from the accoloration vector 
field in the neighborhood of the surface (reference l). 

Sinco the acceleration vector — -, according to the Euler 

dt 

equation, is equal to - — ?rad p and the lattor erpres- 

P 

sion for the case of homogeneity of the medium, with com- 
pressibility considered, can "be rritton equal to - srad 



d-o 

~, the acceleration vector possesses a potontial qp , 



for which there is obtained, "by an integration of the Euler 
equation: 

q> = f(t ) - / ¥ (1) 



Tor the stationary stato, and also for tho nonstat ionary 
case, if the flow at infinity is froe from disturbances , we 
havo in addition, f(tj) = constant. Thus there exists a 
very simple rolation "between the pree-sure field and the 
"acceleration potential." Since the pressure is discontin- 
uous only at the lifting surface and is continuous every- 
where else, the same must hold true for the acceleration 
potential. 

II. flhero tho " linearizod" thoory is enployed (dis- 
turhance velocities evorywhero snail as compared vith tho 
flight velocity V), as is customary in the treatment of 
airfoils, we have t 



"TTeber n'euere Arheiten zur Theor.ie der tra*Senden Flache, 
From 'Proceedings of the Fifth International Congress 
of Applied Mechanics, Oam"orid<;e t Mass., Sept. 1938. 
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— = — + V — (2) 
dt at ax 

or, substituting tho velocity potential , and tho accel- 
eration potential <p: 

«Crad qp = ~ srad 4> + 7 ~ (^rad *) 
6t ox 

from which, interchanging differentiation and integration, 

«p = & + 7 (3) 

at ax 

In tho stationary case — = 0. It is thon possiblo to ob- 

a t 

tain 4 from q> by a simple quadrature: 

.x 

'» (4) 



x 

*(x, y, z) | ^ <p(x«, y, z) dx' 



— oo 



But also in the ^oncral case, $ can "bo computed for ^iven 
«p (x, y, x, t) by an integration of (3). It is nocossary 
to integrate tho acceleration potentials inprossod on oach 
fluid particle: 



0 (x, 7. a. t) - i / <p (x». y. z, t - dx» (4a) 

The shapo of tho lifting surface z = z (x, y, t) can "be oh 
tainod from tho vertical (downwash) velocity xr = d$/dz, from 

the condition = 2 for the stationary case, or from 

ax V 

~ + v It = ^ for the nonstationary case by a second quad- 

OX v o* » 

rature of the sano type. ?or the nonstationary case there 
is obtained: 
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The arbitrary function F takes care of an arbitrary lift- 
ing above the x,y plane of the wins; leading ed*je, and 
likewise of an arbitrary vertical notion of the latter. 

The vortex sheet does not appear in the above formula- 
tion hut, of course, exists as is readily seen from (4) or 
(4a)} if consideration is <?iven to the fact that <p is 
discontinuous at the airfoil and hence in the ahove inte- 
gral expressions, these discontinuities also show up "behind 
the airfoil. 

III. In our considerations thus far, the -wn n I innl ' 
has not been taken into account. In the linearized theory, 
it follows from equation (3) that for reasons of continuity 
the same differential equation that must he satisfied hy <b . 
must also hold for cp . For the incompressible medium, we 
thus have, simply: 



For the oonpres sihle medium we have in the stationary case, 
if the ahove linearization is a^ain applied, 



where c is the speed of sound. 

With the aid of (5) or (5a). the solution can now he 
built up by the "doublet distribution" principle known in 
the electrical theory, A source potential 4>^ can first 

be obtained by distributing sources proportional to the in- 
tensity of the desired 'pressure discontinuity at the air- 
foil- or at the neighboring portion of the x,y plane where, 
for each source at a point x», y', 0, .there is to be sub- 
stituted the singular solution corresponding to differen- 
tial equation (5)' or '(5a); thus, in' the incompressible case, 
the solution const/r, where 



The potential of the source distribution on the unit sur- 
face then becomes 



A cp = 0 



(5) 






(6) 
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where Q is the source intensity. 

To obtain the Telocity potential $, it is 1 necessary 
to pass fron the source to the dipole, which can "be done 
by 

*=-^S (7) 
dz 



IV. Unfortunately, in nany applications of the above 
equations, difficulties arise in carrying out the integra- 
tions, so that it is neceisary to proceed "by the converse 
method, nanely , to find "boundaries of the lifting surface, 
for which conputed source potentials are already available. 
An arbitrary source distribution can then bo built up in 
the forn of a series of such source potentials. 

In the incompressible case, tho procedure has actually 
been applied for the airfoil with circular plan forn. The 
very conplete solution was obtained by W. Sinner (reforonce 
2), who applies olliptic coordinates so that T\ is con- 
stant on confocal ellipsoids of revolution, and li is con- 
stant on the corresponding confocal hypcrboloids (see fig, 
l), « denoting the azimuth. In these coordinates, the 
equation Acp = 0 is satisfied for all 

n n, . m . n n 

cp = P (|a) 0 (ill) coc n*. whore P and Q 
n n to. n n 



denote, respectively, tho spherical harmonics (associated 
Le^ondre functions) of tho first and second kind. If n + m 
n 

is odd, cp^ has a discontinuit7 of a typo "useful for our 
purpose in the x, y plane within tho circle. 

With the above functions which tend to zero on the 
boundary of tho circle, and hence «;ive finite velocities 
there, it is possible to treat cases with impact-free en- 
trance of the flow but not, for oxanple, the case of the 
flat, circular disk set at an an^le to the flow direction - 
for which case, infinite velocities arise at the leading 
el?e. For this purpose it is nocossary to employ additional 

special functions qp t , which are obtained fron ^+3. by a 
suitable differentiation process. We have 

/ a \ n/a 
<P_ = fi M -g H = U ) cos n* 
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The cp n functions become infinite oror the entire boundary 

of tho circle and therefore do not satisfy the flow condi- 
tion at the trailing edge where, according to this condi- 
tion, the disturbance velocity must remain finite. It is 
possible , howevor, froc the infinitely . nany functions 
m 

<P n and <p n to find such a linear combination that, first, 

the vertical velocity w is constant in the interior of 
the circle and, secondly, that the contributions at tho 
trailing edge which become infinite, balance each other, as 
is possible with the Fourier series. By a combination of 
the solution for the flat circular disk with the previous 
solutions, the anglc-of-attack variation of any cambered 
circular surface nay also be treated. 

The results of this theory have been checked in the 
wind tunnel for lift, drag, and pressure distribution by 
M. Hansen, for flat circular plates, spherical segments, and 
an S-cambered surface with fixed center of pressure. The 
theoretical drag tfr is, according to the Hunk stagger the- 
orem, identical with that of the loaded line. The integral 
of the pressure distribution over the loaded surface, how- 
ever, gives a greater drag W s V7hich, on adding the suc- 
tion force at the leading edge, reduces to the value of TT^. 
Since at the ed^es of thin plates the suction force cannot 
actually be fully developed, the value of the true drag 
lies between W L and ¥3 , lying nearer the one or the 
other value, according to the degree of rounding of the 
edge. The surface with fixed center of pressure (fig. 2) 
has been so designed by Kinner that, at some angle of at- 
tack the flow at the entire leading edge is without impact. 
At this angle of attack, therefore, ff s and W L agree. 

•The entire theory holds only for vanishingly small an- 
gles., of attack. Deviations of increasing; magnitude are 
therefore to be expected with increasing angle of attack. 
These appear clearly in the pressure distribution and,' nat- 
urally, also in the values of the lift and drag. If all 
the circumstances mentioned are taken into account and con- 
sideration is also given to the fact that actually the test 
surfaces can only bo designed with finite thickness, where- 
as the theory assumes infinitely thin airfoils, the agree- 
ment of the results given in figures 3 to 5 may be regarded 
as satisfactory, 

V. An investigation in which the problem for the air- 
foil with elliptic plan form is to be solved with the aid 



6 



NACA Technical Memorandum No. 962 



of ellipsoidal .harmonic functions, is. now being conducted 
and likewise, a. computation on the nonatat ionary problem of 
.the vertically flapping circular disk. In the first prob- 
lem, the lack of certain tables of functions led to diffi- 
culties, while in the second tho computation appears to 
proceed smoothly, at least to a first approximation for 
small frequencies. 

The stationary solution for the rectangular airfoil 
is naturally also of importance. This problem has previ- 
ously "been investigated "by H. Blenk (reference 3) "by a dif- 
ferent method which, however, is applicable only to large 
aBpact-ratio wings ( small win? chord), for which ,in most 
applications the loaded-line theory ia found to he -suffi- 
cient. It would therefore ; be of considerable interest to 
obtain a theory for the rectangular airfoil cf small aspect 
ratio,. Unfortunately", however, it appears .that the ex- 
pressions in equation^ ■ (6) and (7) give rise to insurmount- 
able integration difficulties. The author has therefore, 
directed, a computation to be made which is based on the 
principles. of the old airfoil theory and makes use of a 
large number of loaded lines lying one behind the other. 
The case of. four auch loaded lines wasfirst computed, the 
lift distribution of ' each- of these line's being slven by a 
three-term expression and the kinematic condition dz/dx = 
w/V, being satisfied for discrete points between the lines. 
A side investigation showed that the accuracy of. the compu- 
tation, becomes particularly good If the lines are each lo- 
cated at l/4 chord of the surface strips, Into which the 
loaded surface is divided, and the points at which the kin- 
ematic condition Is satisfied are chosen at 3/4 depth of 
the strips. Our coworker, K. TTieghardt, computed the square 
plate as a numorical example, and found that tho distribu- 
tion did not deviate much from the elliptic on any of the 
loaded linos. On tho basis of this result and on tho as- 
sumption that the spanwi.se lift distribution was accurate- 
ly elliptic, ho also treated tho problem for "an Infinite 
number of loaded linos, whore it was now required that .'the 
kinematic condition w = const, bo satisfied in the center 
section only. For tho cho'r'dwisd lift distribution, which is 
represented by the velocity discontinuity 'vCx), an inte- 
gral equation .was now obtained, namOly, 

J 1 — ~ *%=p^ V + - - 

-'■■■>■ 2tt V V sin a 
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for all x f between" +1 and.— l t . A "belli*? tho aspect ratio 

"b/**- 1 t-e . elliptic integral. Tho eolutibn of this in- 
tegral equation was possible only by numerical methods. 
Tor .an expression nade up of four "Birnbaum funfe- 1 

tions" was choson and the integral equation satisfied at 
four points. The required, quadratures had to he carried 
out numerically. This computation save for the various 
aspect ratios the pressure distributions at the center sec- 
tion shown in figure 6. The agreement of the lifts deter- 
mined from these pressure distributions with the test re- 
sults was quite satisfactory, particularly at the small as- 
pect ratios. 

71. With compressibility taken into account and for 
subsonic speeds, the analogy already ^ivon by the author, 
(reforonco 4) nay be directly applied to the preceding re- 
sults obtained for tho lnconpressiblo flow. It is necos- 
sary to write in fornula (6), instead of the provious value 
of r. tho value 

(x - x») + (l - (<y - ;ri) + z ) 



hence the spaco filled with the compressible flow can be 
set to correspond with an affine space filled with a corre- 
sponding incompressible flow "by reducing all dimensions in 
the direction of tho x axis; for example, in the ratio 

1 - 7 /c , while .tho dimensions in the y and z direc- 
tions remain unchanged (Prandtl-G-lauort rule). 

Tho case of supersonic velocity was treatod in detail 
by H. Schlichting (reforonco 5) by the ndw methods. It was 
hore possible to solvo completely the case of tho flat rec- 
tangular plate - with tho restriction, howover, that tho 
disturbance regions which spread out at both side edges un- 
der the Uach angle do not overlap. on the surface. Tho con- 
siderations which have led to formula (6) can here bo ap- 
plied again by sotting in the formula 

r .= /uV x"-»*. - - l) (y - y') 8 + a** 

The values are real only within a doublo cone with tho vor- 
tex at i*, y>,- 0, while outside they are imaginary. Tho 
physical sense requires, howover, that Q(x», y*) should 
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contribute only in the after cone, while in the forward, 
cone and in the outside region l/ r must he set identical- 
ly equal to zero. In order that the forward cone should 
drop out, it is necessary to add a factor 2 to the after 
cone* 

The supersonic problems are simpler than the incom- 
pressible flow and subsonic problems, in that at the lead- 
ing edges no flow arises with infinite volocity (and hence 
no, suction force); also, at the- trailing edges no pressure 
difference approaching zoro is required. Practical partial 
solutions even with the simplest .integral functions are 
therefore obtained for the distribution of the lift density 
Thus Sehlichting, for example, has computed the velocity 
fields for the uniformly loaded rectangular, triangular, 
and trapezoidal wings. The trapezoidal wing constitutes an 
important preliminary work for the theory of the flat rec- 
tangular plate set at angle a to the flow direction. In 
this case the following situation is obtained (fig. 7). 
In a trapezoid ABPE, whose sides AE and BP are in- 
clined at the Mach angle 6, the lift density is constant, 

being equal approximately to 4' ct tan 8 pV 8 /2 ' (whore the 
liach angle a is given by sin a = c/V). Toward the side 
edges the drop takes place in such a manner that along each 
straight line inclined, at an angle cp, for example, BG , 
the lift density is constant. The entire lift distribution 
can therefore be built up by a superposition of a large 
number of uniformly loaded trapezoidal wings with various 
angles cp. In the limit this leads to an integral equation 
which - by a substitution found by I . Lotz - can finally be 
brought to the following form: Lot * = tp.ncp/tan a, and 

f ( be the lift density measured at tho angle cp in the 
t rapezoid-shaped center field. Then 



6=1 




with 



6=o 




The integral equation is to be solved for the boundary con- 
ditions f(4) = 0 for * = 0 and f U) = 1 for * = 1„ 
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The solution has boen' carried out numerically "by Miss I. 
Lot's and J# Pretsch.ln referenco 5. Figure 7 c ;ivos the re- 
sults of the computation. 



SUMMARY 



The "basic idoas of a now nothod for treating tho prob- 
lem of the airfoil are presentod, and a review is *?iven of 
tho problens thus far computed for Incompressible and super 
sonic flows. Test results aro reported for tho airfoil cf 
circular plan fom and tho results are shewn to a^roe well 
with tho theory. As a supplement , a theory based on the 
older methods is prosonted for the rectangular win<; of 
snail aspect ratio. 



Translation by S . RoIbs, 
National Advisory Conmittoe 
for Aeronautics. 
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Figs. 1,6 



M. = 1 




M, = -1 

Figure 1.- Elliptic coordinates 




Figure 6.- Lift distribution in the 

center section of large 
chord flat rectangular plates, after 
K. Wieghardt ( A = aspect ratio). 
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Tigs. 2,3.4,5,7 




Drag obtained 
at c a =0 




for diskj 
HI 

Theoretical j! 
curves of HlnnerT 



Jigura 2.- Surface of circular 

plan form with fixed 
cantor of pressure, after W. Kin- 
ner (flow from the left). 
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Figure 5.- Polar curves of flat 
circular disks with 
rounded leading edge. 
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Figure 3.- Pressure distribution 

over flat circular disk 
set at an angle(the measured values 
are divided "by the dynamic pressure 
q. and tan ot ) . 



Figure 7.- Pressure distribution 
on the flat rectangu- 
lar plate at supersonic speeds, 
after H. Schlichting and I. Lots. 
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Leading edges 




Tigurs 4.- Lift 
coefficients 
of circular 
airfoils 
of various 
camber and 
rounding 
of the 

leading edge 
as a function 
of the angle 
of attack. 
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